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ABSTRACT 

We present a new model-independent method of comparison of near-infrared (NIR) visibility data 
of young stellar objects (YSOs). The method is based on scaling the measured baseline with the 
YSO's distance and luminosity, which removes the dependence of visibility on these two variables. We 
use this method to compare all available NIR visibility data and demonstrate that it distinguishes 
YSOs of luminosity L* < 10 3 i© (low-L) from YSOs of > 1O 3 L (high-L). This confirms earlier 
suggestions, based on fits of image models to the visibility data, for the difference between the NIR 
sizes of these two luminosity groups (Eisner et al. 2004; Monnier et al. 2005). When plotted against 
the "scaled" baseline, the visibility creates the following data clusters: low-L Herbig Ae/Be stars, T 
Tauri stars, and high-L Herbig Be stars. The T Tau cluster is similar to the low-L Herbig Ae/Be 
cluster, which has ^7 times smaller "scaled" baselines than the high-L Herbig Be cluster. We model 
the shape and size of clusters with different image models and find that low-L Herbig stars are the 
best explained by the uniform brightness ring and the halo model, T Tauri stars with the halo model, 
and high-L Herbig stars with the accretion disk model. However, the plausibility of each model is 
not well established. Therefore, we try to build a descriptive model of the circumstellar environment 
consistent with various observed properties of YSOs. We argue that low-L YSOs have optically thick 
disks with the optically thin inner dust sublimation cavity and an optically thin dusty outflow above 
the inner disk regions. High-L YSOs have optically thick accretion disks with high accretion rates 
enabling gas to dominate the NIR emission over dust. Although observations would favor such a 
description of YSOs, the required dust distribution is not supported by our current understanding of 
dust dynamics. 

Subject headings: accretion, accretion disks — circumstellar matter — instrumentation: interferom- 
eters — radiative transfer — stars: formation — stars: pre-main-sequence 



1. INTRODUCTION 

Young stellar objects (YSOs) are surrounded by dust 
and gas leftovers from the process of star formation. Ob- 
servations have revealed that this material concentrates 
into a (protoplanetary) disk accreting toward the central 
star. The disk undergoes the process of dust and gas 
coagulation and formation of larger objects, eventually 
resulting in a planetary system. Hence, YSOs are of a 
topical astronomical interest because they can help us 
understand how planetary systems form and evolve. 

Dust in the disk is very efficient in absorbing the stellar 
radiation. It makes the disk appear optically thick until 
the dust is either removed or coagulated into larger ob- 
jects. Dust is also a very efficient tracer of protoplanetary 
disks in infrared, where dust reemits the energy that it 
absorbed from the star. Direct imaging has proved to be 
especially useful in unraveling properties of YSOs. De- 
pending on the wavelength, images capture different tem- 
perature zones and optical depths of the protoplanetary 
dusty disk. From the terrestrial planet formation per- 
spective, the most interesting is the inner few AU of the 
disk. In this zone dust temperatures can reach a thou- 
sand degrees or more. At these temperatures dust cither 
sublimates or becomes heavily thermally processed, and 
it emits in the near- infrared (NIR). Unfortunately, the 



angular size of this zone is on the miliarcseconds scale 
and until recently unresolvable. 

Thanks to advances in NIR interferometry (e.g. Mon- 
nier 2003), the inner regions of many YSO disks have 
been resolved to date. One of the first surprises coming 
out of these observations is a significantly larger NIR size 
than previously predicted from theoretical models of ac- 
cretion disks (Monnier & Millan-Gabet 2002; Tuthill et 
al. 2001). This led to revisions of the existing theories 
of disk structure by invoking an optically thin inner disk 
hole around the star. The hole is clear of dust because 
of its sublimation, while the gas is unable to provide a 
considerable optical thickness. On the other hand, the 
spectrum shows that the inner disk emits much more 
NIR flux than expected from a simple flat disk. This 
directly implies that the disk structure has to flare up 
in the zone where it reaches its dust sublimation tem- 
perature (T s ^1500-2000K) in order to increase its NIR 
emitting area. Such a model was described by Dullc- 
mond, Dominik, & Natta (2001), where the disk is verti- 
cally puffed up at the inner disk rim because of the direct 
stellar heating of the disk interior. 

Unfortunately, the NIR interferometric observations 
were capable of providing only the characteristic size of 
the emitting region, but they were not good enough to 
constrain the exact geometry of this region. It is only 
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the most recent observations that manage to make a step 
further and detect deviations from simple centrally sym- 
metric images (Monnier ct al. 2006). These deviations 
turned out to be surprisingly small, unlike strong bright- 
ness asymmetries predicted by inclined disk models of 
Dullcmond, Dominik, & Natta (2001). Such a result fa- 
vors a derivative of this model where the inner disk rim 
does not create a vertically flat wall, but rather a curved 
surface (Isella & Natta 2005). In addition, it is still not 
clear from the data if an additional dusty wind compo- 
nent coexists with the disk and contributes to the images 
and NIR flux. According to Vinkovic et al. (2006), such 
an optically thin halo around the inner disk can com- 
pletely explain the NIR flux without a need for the disk 
puffing. 

Another major result coming out of the NIR interfer- 
ometric data is the dependence of inner disk size on lu- 
minosity (Monnier & Millan-Gabet 2002; Monnier et al. 
2005). The size was derived by fitting ad hoc uniform 
brightness ring image models to the measured visibili- 
ties. The inner ring radius Ri n , which is assumed to be 
the dust sublimation radius of the disk, is plotted against 
the luminosity L* of the observed object. The expected 
trend R m cx l}J 2 is followed by YSOs of L+ < 1O 3 L , 
with sublimation temperatures between 1000 K and 1500 
K. Objects with > 10 3 L© deviate from this trend by 
having considerably smaller inner radius than expected. 
Monnier et al. (2005) argue that high luminosity objects 
have optically thick gas within the inner dust-free disk 
hole. This would shield the dust from direct stellar heat- 
ing and enable the dust to survive closer to the star, 
resulting with a smaller inner ring radius. 

The basic approach in these studies is to use predefined 
theoretical models and correlate the object's luminosity 
with a model parameter derived from fitting the model to 
the visibility data. This inevitably leads to questions on 
the validity of the model and conclusions derived from it. 
This is especially true for models of inner disk, consider- 
ing large uncertainties about the inner disk geometry and 
dust properties. In this paper we propose an approach 
that can detect visibility dependence on luminosity with- 
out invoking any model of dust geometry. The approach 
is based on scaling the baseline value such that all ob- 
jects seemingly appear located at the same distance from 
us and have the same luminosity. 

We collected all available NIR visibility data on YSOs 
(13 T Tau, 27 Herbig Ae/Be and 4 FU Ori stars) in the 
literature and in §2 we compared their scaled visibilities. 
We confirm that the circumstellar environment of 
objects with L± > 10 3 i Q differs from less luminous 
objects, with high luminous objects showing smaller 
structures. After that in §3 we plot scaled theoretical 
visibilities of uniform brightness ring, dusty halo, and 
accretion disk over the scaled measured visibilities and 
explore the range of model parameters that are needed 
for accommodating all the data. In §4 we present an 
extensive discussion of various observational aspects of 
YSOs, with a special attention given to the differences 
between low- and high-luminous YSOs. We try to 
sketch a self-consistent model that would explain various 
observed properties of YSOs. The summary of the paper 
is given in §5. 



2. YSO VISIBILITIES WITH SCALED BASELINES 

Studies of the relation between YSO luminosities and 
their NIR visibilities have been based on fitting the con- 
stant surface brightness ring image model to the visibil- 
ity data. The obtained inner ring radius is then corre- 
lated with the luminosity, where the basic assumption is 
that this radius is equivalent to the dust sublimation dis- 
tance from the star (Monnier & Millan-Gabet 2002; Eis- 
ner et al. 2004; Monnier et al. 2005; Akeson et al. 2005a). 
A slightly different approach is used by Vinkovic et al. 
(2006) who use a dusty halo model, but the approach is 
still based on obtaining the sublimation radius from the 
fit to visibility data. Using this approach, Eisner et al. 
(2004) and Monnier et al. (2005) showed that YSOs with 
luminosity >1O 3 L0 have a significantly smaller inner ring 
radius than less luminous objects. 

It is, however, not clear if the ring model is a plausi- 
ble model for fitting the visibility. The NIR disk surface 
brightness is not a uniform ring, but a complicated func- 
tion that critically depends on dust emission and scat- 
tering properties and the disk optical depth structure. 
Moreover, Vinkovic et al. (2006) used a halo model and 
obtained a significantly smaller scatter of inner halo ra- 
dius with luminosity. Here we show that, when explor- 
ing correlations between luminosity and visibility data, 
all uncertainties due to different models of circumstellar 
dust properties and geometry can be avoided by model 
independent scaling of visibility. 

The visibility V(B) is a function of the spatial fre- 
quency B, also called baseline. It consists of two indi- 
vidual additive visibilities: stellar and diffuse. The stel- 
lar contribution to the NIR visibility at currently used 
baselines is a constant because the star remains unre- 
solved. The NIR diffuse part is a contribution from the 
circumstellar dust emission and scattering, while the gas 
brightness can be neglected (we will see later on that this 
is not true for luminous YSOs). 

The diffuse visibility critically depends on two param- 
eters: i) distance to the object, which affects the angu- 
lar size of the image on the sky, and ii) physical size of 
the object, which is controlled by dust radiative transfer. 
The equations of dust radiative transfer have powerful 
intrinsic scaling properties according to which luminos- 
ity and linear dimensions are irrelevant for solving the 
equations (Ivezic & Elitzur 1997). Overall luminosity is 
never an input parameter: only the source spectral shape 
is important. Also, absolute scales of dust densities and 
distances are irrelevant: only the optical depth, geomet- 
rical angles, relative thicknesses and aspect ratios enter 
the equations. Luminosity is important only when we 
need to translate dimensionless radiative transfer solu- 
tions to physical units, where physical dimensions scale 
with yfLl (equation 27 in Ivezic & Elitzur 1997). 

When comparing intrinsic properties of the circumstel- 
lar dust geometry of two objects in the sky, we want to 
make sure that detected size differences are not due to: 
a) differences in distances, and b) size scaling due to 
different luminosities. Therefore, we first need to scale 
objects to the same luminosity and to the same distance. 
Visibility of an image is inversely proportional to the 
characteristic angular size of the image (Ivezic & Elitzur 
1996). In turn, the angular size is inversely proportional 
to distance D, but, as we argue from the radiative trans- 
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Fig. 1. — H-band visibility data with scaled baseline (see equation 1). Letters mark the averaged visibility data from Table 1. The vertical 
error bar is a scatter of measured visibility values, while the horizontal error bar derives from uncertainties in distance and luminosity. The 
upper panel shows the data for a T Tauri star, Herbig Ae/Be stars of L+ < 1O 3 L0 and FU Orionis. The lower panel shows Herbig Be stars 
of L* > 10 Lq. Notice the tendency of low luminosity stars (except for FU Ori) to cluster very compactly, which is a signature of very 
similar images of diffuse radiation, while high luminosity stars have relatively small circumstellar structures. For comparison, dots in the 
lower panel show the location of low luminous stars. Stellar luminosities in units of Lq are indicated together with the stellar name. 



fer scaling properties, proportional to \fhl. Hence, we 
define the scaled baseline 

-Bscalcd = B ■ ^ f —^-i (1) 

D/pc 

which has the same physical units like the ordinary base- 
line, but it removes the intrinsic dependence on distance 
and radiative transfer scaling. 

Although this scaling can be postulated for any system, 
it is meaningful only for systems with a central heating 
source. Binary systems, for example, can have entirely 
different scaling properties. Notice, however, that the 
method is useful for disks with significant accretion lumi- 
nosities because the most of accretion luminosity comes 



from the accretion shock. Hence, the dust in such disks 
is exposed primarily to the sum of stellar and accretion 
luminosity L* + L acc and this sum should replace the 
luminosity factor in equation 1. 

We collected all objects with available H- and K-band 
visibility data and list them in Table 1: 13 objects clas- 
sified as T Tauri stars, 27 Herbig Ae/Be stars and 4 
FU Orionis objects. We excluded from the list objects 
that are most probably B[e] supergiants (MWC 349 A, 
MWC342, HD45677, HD58647) and show them sepa- 
rately only for illustration. Since the goal is to compare 
overall trends in visibility, we average the visibility data 
within bins of about 0.5MA in objects that have many 
measurements taken at similar baselines (the lower and 
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Fig. 2. — Same as Figure 1, but for K-band data. T Tauri, Herbig Ae/Be and FU Orionis are now separated for clarity into their own 
panels. The clustering of low luminosity objects and their distinction from high luminosity objects is clearly visible. The clustering is even 
more evident as we parametrize the data by various visibility functions (see Figures 4-11). 
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upper visibility limits are shown in the table when ap- 
plicable). For a more detailed view of the visibility data 
in individual objects, we suggest following the references 
given in the table. 

Figures 1 and 2 show the visibility data when plot- 
ted against the scaled baseline. The difference between 
low luminosity (L+ < 10 3 L©, low-L) and high luminosity 
(L* > 10 3 L Q , high-L) Hebig Ae/Be objects is now clearly 
visible. It is important to emphasize that this is a model 
independent detection of intrinsic differences between the 
circumstellar environments of low-L and high-L objects. 
Clustering is evident in low-L Herbig Ae/Be stars and 
T Tau stars, while high-L objects cluster at significantly 
larger scaled baselines (except for LkHa 101). The only 
exception are FU Ori objects, which appear more similar 
to high-L Herbig Be stars than to their older counterparts 
- T Tau stars. This suggests possible similarities in cir- 
cumstellar geometry between these two types of objects. 
We address this issue in §3.3.2. 

Figure 2 indicates a slight difference between zones of T 
Tau clustering and low-L Herbig Ae/Be clustering. This 
is a result of differences in the stellar spectral shape, 
where T Tau stars contribute relatively more to the NIR 
flux then Herbig Ae/Be stars. How that influences the 
scaled visibility curves will be more apparent in the next 
section when we plot theoretical models over the data. 

Trends in Figures 1 and 2 are so strong that we can 
use them to identify two anomalous individual cases: 

LkHa 101: Unlike other high-L objects, LkHa 101 
shows a much smaller visibility with the scaled baseline, 
smaller even than low-L objects. This indicates a larger 
NIR emitting area than in other YSOs. This object was 
one of the first YSOs imaged with the NIR interferome- 
try (Tuthill et al. 2001) and it played an important role 
in establishing the existence of the inner disk clearing 
(Tuthill et al. 2002). Prior to the advent of NIR interfer- 
ometry, the canonical model of the inner protoplanetary 
disk was a power-law accretion disk model. Images of 
LkHa 101 reconstructed from aperture masking interfer- 
ometry data showed clearly that this object has a central 
clearing in the disk (Tuthill et al. 2002). The clearing was 
attributed to dust sublimation and LkHa 101 became a 
prototype example of this new disk concept. Our scaled 
visibilities show that this object differs from other YSOs 
studied to date, although the existence of inner clearing 
has been confirmed in other objects. It is possible that 
this object is indeed more similar to the low-L objects 
than high-L objects. There is a large uncertainty in the 
luminosity, distance and evolutionary status of this ob- 
ject. Notice, however, that our scaled baseline in equa- 
tion 1 is not affected by luminosity changes caused by 
changes in measured distance. There has to be an in- 
trinsic error in luminosity to change the scaled baseline. 
For the most recent review of this object see Herbig et 
al. (2004). 

MWC 614'- This objects shows a flat visibility, which 
means it is completely resolved at all used baselines. This 
indicates a much larger structure than other low-L ob- 
jects. Monnier et al. (2006) suggest the possibility of a 
companion star at about 1" distance, which would lead 
to an overestimate of the size of the emission region. 

We excluded from this analysis objects that are most 
probably B[e] supergiants (MWC 349A: Danchi et al. 
(2001); Hofmann et al. (2002); MWC 342: Mirosh- 



nichenko & Corporon (1999); HD45677: Monnier et al. 
(2006); dc Winter & van den Anckcr (1997); HD58647: 
Monnier et al. (2005); Manoj et al. (2002)). Their basic 
stellar parameters are shown in Table 2 and their scaled 
visibilities in Figure 3. When compared with YSOs, 
these objects do not fit into the clustering scheme recog- 
nized in Figures 1 and 2. 

3. THEORETICAL INTERPRETATIONS 

Clustering of scaled visibility indicates very similar im- 
ages, hence presumably similar dust geometry. Varia- 
tions in the model parameters should explain the cluster 
location, size and shape. The advantage of modeling a 
cluster instead of individual objects is that here we can 
address the common properties of these types of objects 
without being misled by potential peculiarities of indi- 
vidual objects. We apply several different types of image 
models and discuss their constrains on the dust geometry 
based on fits to the data clusters. 

3.1. Uniform brightness ring 

A uniform brightness ring emerged as a prototype im- 
age model of NIR interferometry data. It assumes emis- 
sion from a constant temperature disk with a central 
hole, which gives a ring of constant surface brightness. 
The inner radius of the ring is attributed to the dust sub- 
limation, while the outer radius is derived from the to- 
tal ring area required by the measured NIR photometric 
flux. Physical justification for this image configuration 
was found in the puffed up inner disk model (Dullcmond, 
Dominik, & Natta 2001). The most recent NIR interfer- 
ometry measurements of the image departures from cen- 
trosymmetry (Monnier et al. 2006) do not support the 
original version of the puffed up model, but prefer its 
more recent derivative by Isella & Natta (2005). 

The visibility function for a uniform brightness ring 
(along its major or minor axis when inclined) combined 
with an unresolved star is (see Millan-Gabet et al. (2001) 
or Eisner et al. (2004) for details) 



vr s (B) = .n+2 



3^e 2 B{i)) 

tt6 2 B(i) 



Jl(7Tglg(i)) 



(2) 



{B major image axis 

(3) 
B cos(i) minor image axis 

where Ji is the Bessel function and the variables are: 
= fractional contribution of the stellar component to the 
total observed flux at the given NIR wavelength A, i = 
ring inclination angle, and 61,62— inner and outer ring 
size. 

The inner ring radius 6\ is controlled by dust subli- 
mation and grain size. Cooling of smaller grains is less 
efficient than bigger grains, hence bigger grains survive 
closer to the star. Vinkovic (2006) has shown analytically 
that in optically thick disks made of a mixture of grain 
sizes, the inner radius is dictated by the largest grains in 
the mix. They provide shielding of smaller grains from 
direct stellar radiation, which allows smaller grains to 
move closer to the star. The inner disk radius of optically 
thick disks is, therefore, controlled by the largest grains, 
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TABLE 1 

Objects properties and averaged visibility data 
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Distance (pc) 


Luminosity (Lq) 


Baseline (MA) 


H band 


K band 
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T Tau OBJECTS 








LkCa 15 


142±14 (1) 


0.74 (2) 


35.5 




0.95 


2 


BP Tau 


1 42+14 f1 


0.83 (2) 


37.5 




u.yj 


2 


DR Tau 


142±14 (1) 


1.28±0.50 (51) 


38.0 




84-1 03 


3 






48.9 




n &7 n 
u.o 1 -u.yo 




GM Aur 


142±14 (1) 
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2 
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lilt* 1 £j\ 1 \tj J 
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n fifi n Q7 
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u. 1 y u . 00 
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1.4±0.5 (6) 


25.5 
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g 




26.9 




0.97 




RW Aur 


142±14 (1) 


1.7 (2) 


39.0 




0.78 


2 


AS 207A 


165±20 (5) 


2.9±0.2 (4) 


25.5 




0.79-0.96 


4 


V2508 Oph 


165±20 (5) 


4.1±0.8 (4) 


27.0 




0.77-1.00 


4 


DG Tau 


142±14 (1) 


3.62 (2) 


38.8 




n ^7 n fi^i 
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SU Aur 
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U. 1 U-U.oO 












15.0 


fl fi9 Rf> 
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34.8 




4V0 AH 




AB Aur 
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51±14 (12) 


7.5 
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u.yu-u.yo 
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15.0 
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17.5 


n 7o n 77 
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20.0 


0.67-0.77 












22.5 


0.58-0.75 












35.0 




0.41 










39.0 




0.3 




HD344361 (WW Vul) 


550 (21) 


65±5 (18) 


38.2 




0.71-0.84 


6 


V1295 Aql (MWC325, 


290° (12) 


73±10 (10,12,22) 


10.0 


0.95-1.02 




16,13,10,23 


HD 190073) 






14.0 


0.87-0.92 












21.0 


0.68-0.89 












37.0 




0.42 










51.0 




0.00-0.45 




T Ori 


460 (18) 


83 (18) 


38.0 




0.71-0.87 


13 


VV Ser 


350±100 (18,24,25,26) 


43±20 (18,27) 


38.0 




0.54-0.71 


13,19 




45.0 




0.32-0.55 
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Visibility Averages 




Object Names 


Distance (pc) 


Luminosity (Lq) 


Baseline (MA) 


H band 


K band 


Reference (s) 


V1578Cyg (AS477) 


900 (28) 


154±20 (18) 


38.1 




0.79-0.85 


6 


V1977Cyg (Ab442) 


700 (29) 


300±210 (o) 


38.3 




0.30-0.95 


6,13,19 








49.0 




0.45-0.90 




MWC419 (V594 Cas) 


650 (18) 


330 (18) 


39.5 




0.60-0.66 


23 


MWC614 (HD179218) 


255±55 (12) 


\m±Zh h (10) 


11.0 


0.65 




16,23,10 








14.0 


0.62 












22.0 


0.62 












40.0 




0.0-0.45 










51.0 




0.0-0.45 




MWC120 (HD37806) 


360±130 (12,30) 


100±68 (12,30) 


38.0 




0.35 


13 


LkHo 101 


340 c (31,32) 


9000±5700 c (31,32) 


0.5 


0.95-1.05 


0.95-1.0 


32 




1.0 


0.85-0.95 


0.85-0.95 










1.5 


0.75-0.85 


0.65-0.75 










2.0 


0.60-0.70 


0.45-0.55 










2.0 


0.45-0.55 


0.30-0.40 










3.0 


0.30-0.40 


0.15-0.25 










3.5 


0.20-0.30 


0.05-0.15 










4.0 


0.15-0.25 






V380 On 


A A 1 1 r /"IO n 

445±15 (18,30) 


101±16 (18,33) 


11.0 


0.96 




16 




24.5 


0.82-0.88 






MW0147 (HD259431) 


o r\c\ i "i o \ 

800 (18) 


6300 (18) 


21.0 




0.90 


16,8,23 








22.0 


0.98 


1.00 










23.5 


1.00 












39.5 




0.75-0.78 










50.0 




0.71-0.76 




V1685Cyg (MWC340, 


1000 (18,19,34,35,36) 


5000±2000 (18,37) 


23.0 


0.75-0.85 




6,13,16,19 


BD+40 4124) 






38.0 




0.59-0.74 










49.0 




0.39-0.84 




MWC297 (NZ Ser) 


250±50 (38) 


33000 d ±13000 (10,18,38) 


7.5 


0.86-0.92 




10,13,16,23,39 








10.0 


0.81-0.97 












12.5 


0.71-0.81 












15.0 


0.67-0.70 












17.5 


0.63-0.74 












20.0 


0.50-0.67 












21.0 




0.39-0.58 










39.0 




0.00-0.45 










51.0 




0.00-0.45 




MWC1080 


1600±600 (18,40) 


104000±76000 (10,18,27) 


5.0 


0.95-1.00 




10,13,16,19 




7.0 


0.97-1.09 












16.0 




0.93-0.96 










17.0 


0.75-0.88 












20.0 


0.72-0.84 












38.0 




0.45-0.55 




MWC166 (HD53367) 


1150 (18) 


145000±95000 (18,10,41) 


10.0 


0.84-1.00 




10,16 








12.0 


0.87-0.92 












17.0 


0.80-0.92 






Z CMa A 


1100±50 (11,43,44) 


333000±268000 (42) 


27.4 




0.42 


6 




29.0 




0.43 








FU Ori OBJECTS 


2 








FU Ori 


450 (44) 


420±80 (44,46) 


9.6 




0.88-1.01 


47,48 








12.3 


0.93-1.06 












15.9 




0.94-1.03 










21.3 


0.88-1.03 












38.6 




0.86-0.95 










41.5 




0.91-0.96 










47.5 




0.80-0.89 










52.2 


0.86-0.92 












63.2 


0.89-0.93 






V1515 Cyg 


1000±200 (44) 


175±75 (44,46) 


36.7 




0.84-0.92 


49 


V1057 Cyg 


550±100 (50) 


525±275 (44,46) 


38.5 




0.77-0.89 


23,49 


ZCMa SE 


1100±50 (11,43,44) 


510±90 (44,46) 


26.3 




0.40-0.42 


49 
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TABLE 1 
Continued... 



Note. — References: (1) Wichmann et al. (1998), (2) Akeson et al. 
(2005a), (3) Akeson et al. (2005b), (4) Eisner et al. (2005), (5) Chini 
(1981), (6) Monnier et al. (2005), (7) Colavita et al. (2003), (8) Akeson 
et al. (2000), (9) Akeson et al. (2002), (10) Monnier et al. (2006), (11) 
Herbst et al. (1982), (12) van den Ancker et al. (1998), (13) Eisner et al. 
(2004), (14) Meeus et al. (2001), (15) Perez et al. (2004), (16) Millan- 
Gabet et al. (2001), (17) Hernandez et al. (2004), (18) Hillenbrand et al. 
(1992), (19) Eisner et al. (2003), (20) Millan-Gabet et al. (1999), (21) 
Friedemann et al. (1993), (22) Ackc et al. (2005), (23) Wilkin & Ake- 
son (2003), (24) Chavarria-K. et al. (1988), (25) dc Lara et al. (1991), 
(26) Strom et al. (1974), (27) Ackc & van den Ancker (2004), (28) Lada 
(1985), (29) Tcrranegra ct al. (1994), (30) Warren & Hesser (1978), (31) 
Herbig ct al. (2004), (32) Tuthill ct al. (2002), (33) Testi ct al. (1998), 
(34) Lorenzetti ct al. (2002), (35) Shevchenko et al. (1991), (36) Strom et 
al. ( 1972) , (37) van den Ancker et al. (2000) , (38) Drew ct al. (1997) , (39) 
Malbct et al. (2006), (40) Levrcault (1988), (41) Bcrrilli ct al. (1992), 
(42) van den Ancker et al. (2004), (43) Claria (1974), (44) Hartmann & 
Kenyon (1996), (45) Mannings & Sargent (1997), (46) Sandell & Wein- 
traub (2001), (47) Malbct ct al. (2005), (48) Malbet et al. (1998), (49) 
Millan-Gabet et al. (2006b), (50) Straizys ct al. (1989), (51) Muzerolle 
et al. (2003), (52) White & Ghez (2001) 

a Only minimum values are known. 

''van den Ancker et al. (1998) suggest 316^q3-Lq, but we adopt the 
most recent estimate of 1OO±35Zj0 by Monnier et al. (2006). 

c Hcrbig et al. (2004) suggest 700pc, but wc use ~340pc by Tuthill 
et al. (2002) who derived this distance from the analysis of the binary 
companion's apparent motion (luminosity scales with the distance). 

^Significant differences in luminosity reported by different authors. We 
adopt values from Drew et al. (1997). 

e The listed luminosities of FU Ori objects are dominated by accretion 
luminosity. 



TABLE 2 
B[e] STARS 



Object 


Distance 


Luminosity 


Visibility 




(pc) 


(£©) 


Reference 


HD58647 


295±65 (1) 


295±50 (2) 


2 


HD45677 


1000±500 (3) 


14000±7000 (4) 


4 


MWC349 A 


1200 (5) 


55000±25000 (5,6) 


6J 


MWC342 


1000 (8) 


31500±16500 (4,8) 


4 



Note. — References: (1) van den Ancker et al. (1998), (2) Mon- 
nier et al. (2005), (3) dc Winter & van den Ancker (1997), (4) Mon- 
nier et al. (2006), (5) Cohen ct al. (1985), (6) Danchi ct al. (2001), 
(7) Hofmann et al. (2002), (8) Miroshnichenko & Corporon (1999) 



that is, by their sublimation. The most efficient cooling 
in NIR is archived by big (micron size or larger) grains 
because of their gray opacity in NIR and shorter wave- 
lengths. For that reason, big grains reach the minimal 
possible dust distance from the star. Mid-infrared spec- 
troscopy of YSOs shows that the presence of big grains 
is a typical feature of circumstellar dust in these objects 
(e.g. van Boekel et al. 2003; Przygodda et al. 2003; van 
Boekel et al. 2005). Thus, the inner disk radius in YSOs 
is uniquely defined by gray dust grains. 

On the other hand, dust sublimation does not provide 
a unique solution to the inner radius of protoplanetary 
disks because we lack constrains on the exact density 
structure of these disks. Dust dynamics, growth and 
sublimation can make the most inner part of the disk 
vertically optically thin and allow more efficient cooling 




1 2 3 4 5 

Scaled baseline [MA] 

Fig. 3. — Objects that are most probably B[e] stars with in- 
frared excess and erroneously classified as YSOs in observational 
campaigns. This figure shows how they also differ from YSOs in 
their scaled visibility (see Figures 1 and 2). Images of MWC349 A 
are reconstructed from aperture masking intcrfcrometry and here 
we show the visibility along its major and minor image axes. 



of the disk. This reduces the local diffuse radiation and 
enables the optically thin zone of the inner disk to extend 
much closer to the star than the optically thick part of 
the disk. Hence, the inner disk radius due to gray dust 
of temperature Tn is (Vinkovic 2006) 

R in = 0.0344* (^22£) 2 ^[AU], (4) 

where * depends on the details of the disk structure 
and radiative transfer and can be as low as * <~ 1.2 for 
optically thin inner disks or the maximum of * = 2 in 
the case of entirely optically thick disks. In objects with 
a non-ncgligiblc accretion, should be replaced with 
L± + L acc . The scaled baseline (equation 1) is introduced 
from equation 4 by considering the disk size at lpc and 
1L 

e 1 = 0.0688* (^^j [arcsec]. (5) 

The stellar size 0* follows directly from L*, while the 
outer disk size #2 is derived from the requirement of the 
total flux being a sum of the stellar and ring component 
(i.e., V\(B = 0) = 1) 

el = e\ + ^j-, • (6) 

/^cos(z) Sa(Td) 

where 5$a is the Planck function and is the stellar 
temperature. Now we can explore the range of parame- 
ter values constrained by the visibility clusters. We do 
not discuss this model for high-L objects because no sat- 
isfactory fit is possible. 

3.1.1. Inclination 

The clusters are expected to spread in size because of 
randomization of observed disk inclinations. With equa- 
tion 2 we address only the major and minor axes, which 
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Herbig Ae/Be stars L,<1 3 !^ 



0.8 - 



- 



0.4 - 



0.2 - 



f*=0.4 _ 

f*=0.4 

f'=0.4 

J I I I I I I I I L 



i=60°_ f*=0.2 

i=45° f*=0.2 i=60° 

i= 0° f*=0.2 i= 0° , 





T Tau stars 

f*=0.5 

f*=0.5 



f*=0.2 
f*=0.2 




-I 1 1 1 1 1 1 1 1 1 1 1~ 



12 3 

Scaled baseline [MA] 



12 3 

Scaled baseline [MA] 



12 3 

Scaled baseline [MA] 



Fig. 4. — Visibility models of inclined uniform brightness ring. The data is the same as in Figures 1 and 2. The models assume the ring 
temperature Tjj =1500K and ^ = 2 in equation 5 and T*=10000K for Herbig Ac/Be stars and T*=5000K for T Tauri stars. The fraction 
of the stellar flux /* is set to typical values: 0.4 in H- and 0.2 in K-band for low-L Herbig Ae/Be stars and 0.5 in K-band for T Tauri stars. 
While these values reproduce the shape and spread of the low-L clusters, T Tau stars appear larger than the model. An unrealistically low 
value of /* = 0.2 is needed for T Tau stars. 



are the image extremes, while the orientation of baselines 
can exist in between these two position angles. Nonethe- 
less, observations are usually performed within a range 
of position angles, which increases the chance of covering 
the image extremes. These scans of position angles are 
usually incorporated in the visibility vertical error bars. 

We apply equation 2 to the low-L Herbig Ae/Be cluster 
and T Tauri cluster. We use the dust sublimation tem- 
perature of 1500K and the optically thick disk of \& = 2. 
The stellar temperatures are 10000K for Herbig stars and 
5000K for T Tau stars. The fraction of stellar flux /* at 
different wavelengths depends on spectral shapes of stel- 
lar and diffuse radiation. Instead of letting this be a free 
parameter, we use observed spectral energy distributions 
(SEDs) of YSOs to constrain its value. This way we ac- 
tually fit the observed diffuse flux level (relative to the 
stellar flux level), which is needed for realistic modeling 
of both the circumstellar geometry and the SEDs. Ob- 
servations show that the fraction of stellar flux in low-L 
Herbig Ae/Be stars is in the range /* ~0.2-0.6 in H- 
band and /* -0.1-0.3 in K-band (Millan-Gabet et al. 
2001; Vinkovic et al. 2006), while for classical T Tauri 
stars it is /* ~0. 3-0.7 in K-band (Cieza et al. 2005). 
Therefore, we choose medium values of 0.4, 0.2, and 0.5, 
respectively. 

Comparison of models with the data is shown in Fig- 



ure 4. Herbig Ac/Be clusters are nicely reproduced with 
this model for the range of inclination angles between 0° 
and 60°. The model fails in the T Tau cluster. These 
objects appear larger than expected from the model and 
only unrealistically small (i.e., inconsistent with obser- 
vations) values of f*=0. 1-0.2 can reproduce the cluster 
size. This is consistent with the findings of Akeson et al. 
(2005a), who performed a similar study and concluded 
that some T Tau stars appear bigger than the expected 
dust sublimation radius. There is also no positive corre- 
lation between size and accretion, just on the contrary, 
the opposite trend (decreasing size with increasing con- 
tribution of accretion luminosity) has been reported by 
Akeson et al. (2005a). Hence, we conclude that this 
model cannot provide a satisfactory explanation of the 
T Tau data cluster, and we continue with only address- 
ing model parameters in low-L Herbig Ae/Be stars. 

Interestingly enough, the inclination angles up to 60° 
are sufficient to explain the low-L cluster size even though 
we have not yet varied other parameters. Since we ex- 
pect other parameters to spread the model further, this 
might be an indication of images appearing as if they 
had inclination angles in a smaller range than expected 
from random orientations. We address this issue further 
below. 
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T D =2000, f*=0.6 - - 

T D =2000, f*=0.3 — 

T D =1500, f*=0.6 — 
T„=1500, f*=0.3 




T,= 12000, f =0.6 i=60° f*=0.3, T,= 12000 

T,= 7500, f*=0.3 i=60° f*-0.1 . T,=7500 

T,= 12000, f*=0.6— — — i=0°— — — f-0.3, T,= 12000 

T*= 7500, f=0.3 i= 0° f-0.1, T.=7500 



1 2 

Scaled baseline 



1 2 
Scaled baseline [MX] 



Fig. 5. — Low-L Herbig Ae/Be stars modeled with inclined uni- 
form brightness ring. The model parameters are the same as in 
Figure 4, except for the stellar fraction /* and ring temperature 
T]j. The stellar flux fraction covers the lower and upper limit, 
combined with the lower and upper inclination angle limits from 
Figure 4. Data clusters appear smaller than the range of these 
parameters and with a clear preference for dust temperatures of 
1500K. T Tau K-band cluster cannot be successfully fitted with 
this model (see Figure 4). 



T,= 12000, f - =0.6 i=60° f =0-3. T,= 12000 

T,= 7500, f - =0.3 i=60° f-0.1. T.=7500 

T,= 1 2000, f '=0.6 — — — i=0°— — — f'=0.3, T,= 12000 

T= 7500, f*=0.3 i= 0° f-0.1, T.=7500 

I ... i I i i .. I I i i . . I 




1 2 3 1 2 3 

Scaled baseline [MX] Scaled baseline [MX] 

Fig. 6. — Same as Figure 5, but for different stellar temperatures. 
The disk temperature is set to 1500K. 



3.1.2. Stellar flux fraction and dust sublimation 

Next we explore the above mentioned limits in the stel- 
lar flux fraction. Figure 5 shows models with limiting 
values of /* combined with disk inclinations. Combining 
low /* with low i and high /* with high i brings up ex- 
tremes in visibility functions. Comparison with the data 




12 3 

Scaled baseline [MX] 



1 2 

Scaled baseline [MX] 



Fig. 7. — Same as Figure 5, but now for optically thin disks 
with = 1.2. Although the model of an optically thin disk would 
work for small inclination angles, the data clusters clearly support 
an optically thick disk when a full range of inclination angles is 
considered. 



shows that clusters appear slightly more compact than 
these limits, but this could be just the result of a small 
dataset. However, freedom in choosing the disk tempera- 
ture (i.e., dust sublimation temperature) results in large 
model deviations from the data, as shown in Figure 5. 
The temperature of 2000K is often suggested to explain 
the data in some individual objects, but when it comes to 
explaining the collective dataset, then the temperature 
of 1500K emerges as the most plausible choice, assuming 
that visibility data clusters are indeed a result of very 
similar physical conditions in all these objects. 

3.1.3. Stellar temperature and disk opacity 

Objects comprising the data clusters have a range of 
stellar temperature from 7500K to 12000K for low-L Her- 
big Ae/Be stars. Figure 6 shows how this spread in tem- 
perature affects the model. Since the stellar spectrum 
at these temperatures peaks at wavelengths shorter than 
NIR, changes in the stellar temperature do not have a 
significant effect on the NIR images. Differences of these 
models from the 10000K models in Figure 5 are small 
and conclusions from §3.1.2 also hold here. 

In addition, we explore the possibility of an optically 
thin inner disk and apply \& = 1.2 to these models. Fig- 
ure 7 shows that optically thin disks would be an op- 
tion if disks are almost face on, but not when the whole 
cluster of visibility data is considered. Since \I><~1.2 is 
the minimal possible value, maybe some intermediate 
(1.2 < < 2) values exist. However, the NIR excess in 
Herbig Ae/Be stars cannot be explained by optically thin 
disks alone. Instead, an additional circumstellar struc- 
ture, such as a puffed up disk or a dusty halo, has to be 
invoked to reproduce the observed amounts of NIR flux 
(Vinkovic et al. 2006). 
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Herbig Ae/Be stars 



t v =0.80, T,= 12000 K 

t v =0.10, T,= 12000 K 

T v =0.80, T,= 10000 K 

r v =0.15, T,= 10000 K 

t v =0.80, T,=7500 K 



t v =0.20, T,=7500 K 



T Tauri stars 



t v =1.00, T,=5,000 K 

i\,=0.60, T =5,000 K 



t v =0.20, T,=5,000 K 



0.8 - 



^0.6 - 



> 



0.4 - 



0.2 - 




1 2 3 1 2 3 

Scaled baseline [MA] Scaled baseline [MA] 



12 3 

Scaled baseline [MA] 



Fig. 8. — Visibility models of dusty halos. The models assume a spherical halo for simplicity. The dust is lfim silicate grains sublimating 
at 1500K. Lines show visibility functions for various optical depths and stellar temperatures. Optically thin halos are sufficient to explain 
the data clusters. 



3.2. Dusty halo 

An alternative model to disk geometries is a dusty halo 
around the disk inner regions. Compact (~10 AU) op- 
tically thin halos have been invoked to explain the NIR 
spectrum of Herbig Ac/Be stars and their NIR visibility 
(Johns-Krull & Valcnti 2001; Vinkovic et al. 2006). The 
exact physical mechanism for creating such a halo is not 
known, but the variability of YSOs (discussed in §4.2) 
supports its existence as a dusty outflow. 

In its most rudimentary form, the halo can be approx- 
imated with a spherical geometry. We emphasize that 
this is just a convenient approximation, while a more re- 
alistic description of a dusty outflow probably requires a 
flattened and clumpy halo. Nevertheless, optically thin 
halos are dominated by direct stellar heating and, there- 
fore, will maintain similar temperature profiles, no mat- 
ter the exact geometry. We can also ignore the disk heat- 
ing of the halo because optically thin halos are transpar- 
ent to the disk emission. The exact images of the disk 
and halo model ultimately depend on detailed properties 
of the halo and the intrinsic ratio between the disk and 
halo surface brightness. These images also depend on 
the inclination angle, which is something that we cannot 
address here with our simple approximation. We caution 
that halo models cannot be dismissed by postulating that 
they should produce centro-symmetric images (e.g. Isella 
et al. 2006) . Speckle interferometry images of the Herbig 



Be star R Mon are an illustrative example of compli- 
cated asymmetric images produced by a parabolic dusty 
outflow combined with the inclined disk (Weigelt et al. 
2002). 

The halo radiative transfer is solved with the code 
DUSTY (Ivezic, Nenkova, & Elitzur 1999), which takes 
advantage of the scaling properties of the radiative trans- 
fer problem for dust absorption, emission and scattering 
(Ivezic & Elitzur 1997). The stellar spectra used in our 
modeling are taken from Kurucz models. The silicate 
dust optical constants are from Dorschncr et al. (1995) 
{x = 0.4 olivine). The radial density profile is described 
with r~ 2 since Vinkovic et al. (2006) showed that it can 
reproduce the NIR SEDs of Herbig Ae/Be stars. The 
outer halo radius has to be large enough to extended be- 
yond the distance where the halo brightness drops below 
interferometric or photometric detection. Hence, we fix 
it to be 100 sublimation radii. Visibility functions at 
specified wavelengths are part of the computational out- 
put from the code. The output visibility already has the 
stellar and diffuse components included, thus, we do not 
need to specify /* in this modeling. The variable that 
indirectly controls the amount of diffuse flux in optically 
thin halos is the radial optical depth (see Appendix in 
Vinkovic et al. 2003). 

The provided spatial frequency (fausty from DUSTY is 
the frequency scaled with the angular size of the inner 
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. t v =0.80, 0.2jnm 
t v =0.15, 0.2jnni 

. t v =0.80, Ijitm 
t v =0.15, 1jum 

. t v =0.80, lOjlim 
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12 3 

Scaled baseline [MA] 



12 3 

Scaled baseline [MA] 



Fig. 9. — Same as Figure 8, but now for different grain sizes. The 
data clusters belong to low-L Herbig Ae/Be stars. The model stel- 
lar temperature is 10000K. Lines show variations in optical depth 
and grain size. If dusty outflows exist in these objects, then the 
visibility data suggest dust grain size of ~ l/im. 



cavity (Ivezic & Elitzur 1996). Since the inner cavity 
radius n for = 10 4 L Q is also provided in the output, 
we can derive the scaled baseline from 



Scaled = 1.54 x 10 14 g dusty (—) 1 (—) . (7) 

Comparison between halo visibility models and data 
clusters is presented in Figure 8. Dust grains used in the 
models are 1/nm in size and sublimating at 1500K. We use 
stellar temperatures typical for these stars and find that 
optically thin halos can explain the data clusters. The 
required visual optical depth in low-L Herbig Ae/Be stars 
is Ty ~ 0.15 — 0.8, which is in agreement with ry > 0.2 
derived from fits to SEDs (Vinkovic et al. 2006). The T 
Tau cluster displays a slightly larger upper optical depth 
limit of Ty ~ 1. 

The grain size of ~ lfim seems optimal to fit the data. 
Figure 9 shows how the visibility changes when the grain 
size is reduced to 0.2^m or increased to 10/im. These 
changes are caused by variations in the sublimation ra- 
dius R in of the halo, where larger grains reduce and 
smaller grains increase this radius. Similar grain sizes 
fitting the data mean that the inner halo radii are also 
very similar. 

The model is unsuccessful in explaining the visibility 
of high-L Herbig Be stars. The only way to reach the ob- 
served visibility levels of these stars is to keep the halo at 
tiny optical depths, which creates a flat visibility due to 
an unresolved star at observational baselines. This was 
the approach used by Vinkovic et al. (2006) in fitting 
one high-L star (MWC 297). Since observed visibilities 
show a non-fiat structure, this model is unsatisfactory. 
This shows that if the halo model is correct in low-L ob- 
jects, then high-L objects differ from low-L objects by not 
having a dusty halo, which is not surprising considering 



that large luminosities would impose a large radiation 
pressure on optically thin dusty halos. Such a scenario 
would be in agreement with observations showing that 
YSOs of spectral type earlier than about B5 disperse 
their circumstellar environment much faster than later 
spectral types (Fuente et al. 1998). Hence, the NIR vis- 
ibilities of high-L objects imply that these are stars that 
are still surrounded by circumstellar gas and dust, but 
the process of dispersal has already started within their 
immediate environment. 

3.3. Accretion disks 

Prior to the discovery of inner disk holes in low-L Her- 
big Ae/Be stars produced by dust sublimation (Monnier 
& Millan-Gabet 2002), accretion disks were assumed to 
be the main source of the NIR excess in YSOs. Since ac- 
cretion disks extend closer to the star than the dust sub- 
limation radius, they are still a favored explanation for 
the NIR visibility data of high-L Herbig Be stars (Mal- 
bet et al. 2006; Monnier et al. 2005; Eisner et al. 2004, 
2003; Monnier & Millan-Gabet 2002; Millan-Gabet et al. 
2001). 

3.3.1. Accretion disk model 

The total luminosity of stars with accretion disks is a 
sum of the stellar and the accretion luminosities, L toi — 
L± + L&cc where L acc = GM i ,M &cc /2R i ,. Combining stel- 
lar heating (Friedjung 1985) (for R > 2R+) and viscous 
heating (Lynden-Bell & Pringle 1974) of a geometrically 
thin, optically thick disk results in the disk temperature 
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where 9 is the scaled angular size derived from the disk 
radius R and the total luminosity 
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The inner disk edge size 9 m (radius i?i n ) and disk incli- 
nation angle i are free parameters. The stellar fraction 
of NIR flux f£ can be computed by integrating the disk 
emission 
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where the upper integral limit can be any size where the 
disk temperature drops below detection at NIR wave- 
lengths. 

The visibility is a sum of narrow rings over the disk sur- 
face. Each ring has the visibility described by equation 
2, hence the accretion disk visibility is 
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Fig. 10. — Accretion disk models of high-L Hcrbig Be stars. The data are the same as in Figures 1 and 2. The legend lists parameter 
values used in modeling: stellar temperature T* in K, stellar luminosity and mass M* in solar units, accretion rate M acc in M^/yr, 
and inner disk radius i?= n in solar units. Presented models are for face on disks. Other listed parameters are derived from the model: 
disk temperature T= n in K at the inner disk edge, accretion luminosity L= n in solar units, and stellar flux fraction in K-band (f*(K)) and 
H-band (f*(H)). See §3.3 for more details. 



where AO is a small angular step and B(i) is the scaled 
baseline defined in equation 3. 

Unlike in dust ring models, here we do not have a firm 
theoretical limit on the inner radius i?= n , but if the disk 
accretes onto the star, then the inner radius of several 
stellar radii is expected. The main problem in constrain- 
ing these models is that precise stellar parameters are not 
well known for high-L Hcrbig Be stars. On top of that, 
accretion luminosities are relatively small compared to 
stellar luminosities, which makes accretion rates difficult 
to constrain. 

All these uncertainties are illustrated in Figure 10 
where we compare various accretion disk models with the 
data (we ignore anomalously low visibilities of LkHa 101 
in this discussion, as described in §2). We display mod- 
els with no accretion (a purely reprocessing disk) and 
with a very high accretion rate of 10~ 4 M Q /yr to illus- 
trate how models with almost any accretion rate can be 
build to fit the visibility data (for further discussion see 



e.g. Malbet et al. 2006). Models are calculated only for 
face on inclination, since any loss of flux in inclined disks 
can be compensated by an increased inner disk radius 
(Millan-Gabet et al. 2001). Also, inclination increases 
the visibility, which again calls for a larger inner radius 
to fit the data. Hence, presented models put a lower limit 
on i?; n . 

3.3.2. The nature of disk in high-L Herbig Be stars 

It seems that some models with i?= n ~ 5i?* can fit the 
visibility, which would be in agreement with the mag- 
netic accretion radii implied from T Tau stars (Kenyon 
et al. 1996). However, Figure 10 shows that models with 
>10i?* are also applicable. Moreover, a detailed model- 
ing of MWC 297 yields R in - 17i?* (Malbet et al. 2006), 
which raises the question of what is the physical process 
that truncates the disk at these radii. 

Monnier et al. (2005) suggest that the NIR visibility 
of high-L objects is due to dust in the accretion disk. In 
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Fig. 11. — Same as Figure 10, but for FU Ori stars. The K-band 
data are shown. The stellar temperature, mass and luminosity are 
5000K, IMq and 4Lq, respectively. Lower than expected visibil- 
ities are most probably an indication of over-resolved large dusty 
halos (Millan-Gabet et al. 2006b). 



that case the inner radius is dictated by the dust subli- 
mation temperature. Figure 10 shows that sublimation 
temperatures would require R ln > 20. Such disks would 
be optically thick in their purely gaseous part inside of 
the dusty disk hole. This would allow shielding of the 
dust from direct stellar radiation, so that the dust can 
move closer to the star than the ordinary dust sublima- 
tion distance. The extent of this shielding can be esti- 
mated by scaling the high-L visibility clusters in Figures 
1 and 2 to the location of low-L clusters. We find that all 
the stellar flux shorter than ~ 0.6/im has to be removed 
to make this process work. The obtained scaling factor 
is <~7, which also shows how smaller high-L objects are 
from the dust sublimation size typical of low-L objects. 

Optically thin gas can be optically thick in some molec- 
ular and atomic lines, but it would be transparent at 
other wavelengths and incapable of producing such a 
severe reduction of the stellar flux. The shielding is, 
therefore, plausible only if the gas is completely op- 
tically thick, which also implies high accretion rates 
(> 10- 7 M Q /yr) (Muzerolle et al. 2004; Bell & Lin 1994; 
Hartmann et al. 1993). On the other hand, this cre- 
ates a problem because the side effect is a strong near 
IR radiation that can actually overshine the dusty disk. 
Muzerolle et al. (2004) demonstrate this in the case of a 
puffed up dusty disk, where gas emission from accretion 
disks of M acc > 10 _6 M o /yr would dominate the NIR. 
Since shielded dusty accretion disks do not have such a 
puffing, their near IR dust emission is smaller than from 
puffed up disks and, therefore, competes with smaller 
accretion rates. 

The effect is even more pronounced in NIR visibilities 



where a bright optically thick gas would affect images by 
decreasing the observed inner disk radius. Comparison 
between disk models of different inner radii in Figure 10 
shows that any such decrease would significantly alter 
the visibilities. Thus, if dust is the source of NIR visi- 
bility of high-L objects, then the gas inside of the dust 
sublimation radius has to be optically thick for stellar 
radiation, but not bright in the NIR at the same time. 
Scattering of UV by low-density gas was suggested by 
Monnier & Millan-Gabet (2002) as a possible solution, 
but this cannot work because i) large optical depths are 
needed to block almost all of the scattered stellar flux, 
and it) even if all of the stellar UV flux is removed, it is 
still not enough to make high-L scaled visibilities similar 
to low-L visibilities. This dilemma of how to reconcile 
such properties of the gaseous disk within the dust sub- 
limation zone is actually quite old, and it was already 
recognized in a slightly different context by Hartmann et 
al. (1993). 

It is, therefore, possible that instead of dust, the NIR 
flux of high-L objects originates from gas emission de- 
spite confusions about the inner disk radius. This would 
suggest that these stars have high accretion rates and 
their circumstellar environment is evolutionary similar 
to FU Ori objects (young counterparts of T Tau stars). 
Since both types of objects are of a similar age, they 
probably share similarities in the accretion process. Fast 
pre-main-sequence evolution of massive stars actually 
requires high accretion rates (Palla & Stahlcr 2000). 
Also, observed high mass loss rates imply high accretion 
rates in massive stars (Cesaroni et al. 2006; Shepherd 
& Churchwell 1996), similar to massive winds supported 
by high accretion rates in FU Ori stars (Sandell & Wcin- 
traub 2001; Calvet et al. 1993). The main difference is 
the ability of highly luminous Hcrbig Be stars to rela- 
tively quickly disperse their circumstellar matter (Fuente 
et al. 1998). 

The standard viscous disk model predicts proportion- 
ality between the disk mass and the accretion rate (e.g. 
Calvet et al. 2000). Thus, high accretion rates would 
suggest higher disk masses in high-L than in low-L ob- 
jects. Observations, however, show exactly the opposite 
(Fuente et al. 2003). It may be that the dispersion of 
circumstellar matter around high-L stars starts with disk 
erosion due to disk and stellar winds. This may reduce 
the disk mass and explain the observed decrease in disk 
masses. Their destructive nature may also explain why 
these stars form a data cluster at ^7 times larger scaled 
baselines than low-L YSOs instead of producing a range 
of values due to disk optical depth variations between 
objects. Namely, it may be that only high density and 
high accretion disks can survive so close to these lumi- 
nous stars, which would result in disks with very similar 
NIR signatures. 

High accretion rates may also be implied from the re- 
cently observed correlation between the accretion rate 
and the square of the stellar mass of low-L YSOs (e.g. 
Calvet ct al. 2004; Muzerolle et al. 2005; Garcia Lopez 
et al. 2006; Natta et al. 2006, and references therein). 
If this correlation is not just a selection effect (Clarke 
& Pringlc 2006), then it probably extends to high-L 
objects, too. The required accretion rates would be 
M acc > 10~ 7 M Q /yr, which is exactly what we suggest 
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based on the NIR visibilities. 

There are additional complications in deducing the na- 
ture of disks around Herbig Be stars. MWC 166 is, for 
example, a binary with a mean distance of ~1.7AU be- 
tween components and an orbital eccentricity of e=0.28 
(Pogodin et al. 2006). The binary is surrounded by a 
common gaseous envelope and the existence of the disk 
in this system is a more puzzling problem than our gen- 
eral discussion about Herbig Be stars. We also ignore any 
contribution from the NIR free-free emission. Although 
this is a good approximation for Herbig Be stars, it is 
possible that some objects have a significant amount of 
free-free emission (e.g MWC 297: Porter et al. 1998) 

3.3.3. FU Ori accretion disks 

Accretion disks in FU Ori stars are corroborated with 
NIR visibility measurements, as evident from Figure 11 
where we compare the visibility data of FU Ori stars 
and accretion disk models. Modeling is simplified in 
this case due to similarities in stellar properties of FU 
Ori stars. We use face-on disks and stellar tempera- 
ture, mass and luminosity of 5000K, 1M and 4L Q , re- 
spectively. Classical accretion disks of Rin—Z-^R^ and 
M acc = 10~ 4 — 1O _5 M0 /yr reach the upper observed vis- 
ibility levels (see also Malbet et al. 2005, 1998; Millan- 
Gabet et al. 2006b). There is a problem with anoma- 
lously low visibility of Z CMa-SE, which is more con- 
sistent with T Tau visibilities. Other FU Ori stars also 
show a slightly smaller visibility than expected from pure 
accretion disks, which indicates an additional larger com- 
pletely resolved circumstellar structure. Millan-Gabet et 
al. (2006b) argue that this is due to a large dusty en- 
velope. Such a disk+envelope structure was previously 
predicted from the overall infrared spectral energy dis- 
tribution (Hartmann & Kenyon 1996). This configura- 
tion and derived accretion rates were recently confirmed 
by modeling the infrared spectra of FU Ori stars taken 
with the Spitzer Space Telescope (Green et al. 2006). 
However, Quanz et al. (2006) modeled their mid IR visi- 
bilities of FU Ori without any additional structure such 
as a dusty envelope and concluded that the presence of 
accretion disk is sufficient to explain this object. 

4. DISCUSSION 
4.1. The importance of dust dynamics 

We usually approach this kind of study with the as- 
sumption that circumstellar environments of YSOs share 
enough similarities to be successfully described with one 
universal theoretical model. This is a very strong as- 
sumption considering that the spectra reveal very dy- 
namic gaseous disks around YSOs (e.g. Mora et al. 2004). 
Since gravity and gas drag are two forces that dictate 
dust dynamics within protoplanetary disks, dust should 
also display a very dynamic behavior. The infrared spec- 
tra provide evidence that this is indeed the case: the 
mid IR dust features exhibit clear signatures of small 
grains in the disk surface ( Acke & van den Ancker 2004) , 
which is possible only if dust dynamics persistently re- 
plenishes the surface with these grains (Dullemond & 
Dominik 2005). 

Dust dynamics is, however, ignored when it comes to 
formulating dusty disk models used in calculating syn- 
thetic infrared spectra and images. The disk is modeled 



with a smooth surface that seemingly appears dynam- 
ically passive. This is a good approximation when it 
comes to modeling data taken at one epoch, but multi- 
epoch data reveals that YSOs are far more complicated 
than such simple disk models. Basically all pre-main- 
sequence objects are spectrally and photometrically vari- 
able, differing only in the amplitude and rate of variabil- 
ity. This variability includes also the near and mid IR 
wavelengths where dust emission dominates (e.g. Chen 
& Jura 2003; Grinin 2000; Herbst & Shevchenko 1999; 
Skrutskie et al. 1996; Liu et al. 1996; Hutchinson et al. 
1994; Prusti & Mitskevich 1994). 

While some mid IR variability can be attributed to 
variations in the stellar (or accretion) heating of the 
dust, many of these variabilities are not accompanied by 
changes in the luminosity of the central source. Such 
cases are a clear sign of dust dynamics producing mod- 
ulations in the dust emission and/or stellar obscuration. 
On top of that, Vinkovic (2006) showed analytically that 
the dust sublimation zone of optically thick dusty disks 
cannot be constrained purely by radiative transfer. Sub- 
limation leads to vertical optical thinning of the disk, 
which, combined with dust dynamics, leads to entirely 
new radiative transfer solutions. This result demon- 
strates that inferring a realistic geometry of dust dis- 
tribution in the inner disk region of optically thick pro- 
toplanetary disks requires dust dynamics as well as ra- 
diative transfer. 

Although the NIR visibilities are currently observed 
with a limited u-v coverage, dust dynamics responsible 
for strong variabilities should be detectable in multi- 
epoch visibility observations. In addition, the newly 
employed NIR closure-phase measurements enable inter- 
ferometric detections of dusty disk inhomogeneities, as 
recently demonstrated in the case of AB Aur (Millan- 
Gabet et al. 2006a). Time evolution of such brightness 
asymmetries may significantly improve our understand- 
ing of processes responsible for the infrared variability of 
YSOs. 

This is an important topic because current disk models 
have problems in explaining variabilities caused by tran- 
sient dust obscurations of the central star (UXOR vari- 
ables). This variability is manifested as changes in dust 
extinction, accompanied by the "blueing" effect and in- 
creased polarization at minima (Shakhovskoj et al. 2005; 
Rodgers et al. 2002; Rostopchina et al. 2001; Grinin et 
al. 2001; Grinin 2000; Skrutskie et al. 1996; Grinin et al. 
1994; Hutchinson et al. 1994). The photometric minima 
effects are a result of increased relative contribution of 
scattered starlight by dust in the total spectrum (Natta 
& Whitney 2000). 

Temporal properties of obscuration events affecting vi- 
sual and NIR wavelengths indicate that dust clouds usu- 
ally appear in the inner disk regions at or close to the dust 
sublimation zone. If they belong to inhomogeneities con- 
strained to the disk then disk inclinations have to be large 
(Dullemond et al. 2003). The necessesity for such a cor- 
relation was originally suggested by Natta et al. (1997). 
If we assume the most optimistic inner disk puffing of 
H/R -0.2 (Vinkovic et al. 2006) then the line of sight 
would be affected by the inner disk when inclination an- 
gles are i >80°. This is inconsistent with the imaging 
data, which show disks with smaller inclination angles. 
In particular, the UXOR star CQ Tau (Shakhovskoj et 
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al. 2005) has i = 48°t| derived from the NIR visiblities 
(Eisner et al. 2004), i = 33° ±5 from the mid IR imaging 
(Doucet et al. 2006) and i ~ 60 — 70° from mm-imaging 
(Testi et al. 2003), while the UXOR star VV Scr has 
i = 42° t j derived from the NIR visiblities (Eisner et al. 
2004). 

We do not agree with the interpretation by Isella et 
al. (2006) that their derived disk inclinations (from fits 
to the NIR visibility) of 40-55° for CQ Tau and 50-70° 
for VV Ser are in agreement with the dust obscuration 
model. These inclinations would actually support the 
idea of dust clouds being ejected to >9 (CQ Tau) and 
>4.5 (VV Ser) vertical scale heights above the dusty disk, 
which would dynamically decouple them from the disk 
(i.e., their dynamics would not be controlled by gas-drag 
any more). 

4.2. Are dusty halos a plausible option? 

Dust obscuration events are easier to explain if we al- 
low for a possibility of dust clouds being ejected to mod- 
erate polar angles above the disk mid-plane. Notice that 
this would still imply a dependence of variability on the 
inclination angle, which can explain the observed weak 
correlation between inclinations measured in v sin i and 
the amplitude of photometric variability (Grinin & Ko- 
zlova 2000) . Similar correlation between the observed po- 
larization and v sin i does not exist, but polarimetry sup- 
ports the idea of a flattened halo above the disk (Yudin 
2000). In addition, numerical models reproduce pho- 
tometric and polarimetric variability of UXORs using 
the ad hoc assumption of obscuring dust clouds above 
the scattering dusty disk (Natta & Whitney 2000). We 
note that the variable accretion luminosity model put 
forward by Herbst & Shevchcnko (1999) in order to ex- 
plain the variability of UXOR objects was successfully 
challenged by the proponents of dust obscuration model 
(Rostopchina et al. 2001; Grinin et al. 2001; Grinin 2000). 

Once a cloud is out of the dense gaseous disk, gas drag 
force becomes negligible and radiation pressure takes 
over. The cloud would be eventually blown away, cre- 
ating a dusty halo-like outflow above the disk. Changes 
in such an outflow would manifest itself by changes in the 
near and mid IR thermal emission. This IR variability 
is not dependent on the inclination angle, assuming op- 
tically thin emission. It would also be uncorrelated with 
the stellar spectrum variability. Although simultaneous 
visual and infrared observations are rare, such variabil- 
ity events have been documented (Chen & Jura 2003; 
Eiroa et al. 2002; Prusti & Mitskevich 1994; Hutchinson 
et al. 1994). Particularly interesting is the case of HD 
163296 in which Sitko ct al. (2004) detected an event of 
a major increase in the near and mid IR emission. Since 
HST imaging shows a Herbig-Haro flow from this star 
(Grady et al. 2000a), they conclude that the observed 
infrared "flare" is possibly an ejection of a forthcoming 
Herbig-Haro object. If correct, this would directly link 
stellar outflow with the infrared variability due to dusty 
outflow. 

Unfortunately, despite these seemingly convincing ar- 
guments, the existence of a dusty outflow is not straight- 
forward and obvious. The biggest drawback in halo mod- 
els is the lack of a known force capable of lifting a dust 
cloud out of the disk. The common assumption in cir- 



cumstellar disk modeling is that gravity and gas drag 
do not allow dust to exit the high density gaseous disk. 
Hence, the existing studies of NIR visibilities are domi- 
nated by disk models. Circularly symmetric images are 
sometimes invoked as a feature of dusty outflows not in 
agreement with NIR interferometry observations. How- 
ever, in §3.2 we argue that this is an ill-formed argument 
because halos can certainly display asymmetric images. 

There was great excitement for some time about the 
possibility of cometary activity in Hcrbig Ae/Be stars 
detected in spectral line variations. This is called the (3 
Pictoris phenomenon, after the same effect first observed 
in the dust debris disk of (3 Pic (Grady et al. 2000b). Ac- 
cording to this interpretation, accretion episodes accom- 
panied by redshifted absorption components in a num- 
ber of metal lines are signatures of infalling evaporating 
comets. This would nicely fit into the above scheme of 
transient dust clouds. Unfortunately, further observa- 
tions showed that the composition of infalling gas is con- 
sistent with the accretion disk gas (Natta et al. 2000), 
while numerical models revealed that spectral signatures 
of comets cannot be detectable due to suppression by 
strong stellar winds in Hcrbig Ae/Be stars (Beust et al. 
2001). Since then there have been no other suggestions 
for a process that could lift a dust cloud out of the disk. 

A new imperative to explore dusty halos may come 
from the NIR visibility studies. The large NIR excess 
(often called the NIR "bump") of Herbig Ae/Be stars re- 
quires a vertical disk puffing at the inner disk edge (Isella 
et al. 2006; Dullcmond, Dominik, & Natta 2001). The 
NIR images based on this model have a large skewness - 
the brightness ratio between two opposite sides of the im- 
age centered on the star. Monnicr et al. (2006) recently 
measured the NIR interferometry closure phases that can 
reveal the amount of skewness and found that observa- 
tions do not support model predictions 1 . However, their 
claim is based mostly on high-L objects that do not ex- 
hibit visibilities consistent with the puffed up disks any- 
way (see §3.3.2). The main question is whether low-L 
Herbig Ae/Be stars can be described with this model or 
not. It turns out that closure phases are inconclusive for 
their skewness determination due to insufficient telescope 
resolution. 

On the other hand, detailed modeling of the NIR visi- 
bilities with the puffed up disk model is not successful in 
reproducing strong near IR bumps (MWC 758 and VV 
Ser: Isella et al. 2006). It was exactly these objects with 
a strong NIR bump that Vinkovic et al. (2006) point 
out as the most puzzling and most difficult to explain. 
The model also has problems with reconciling variability 
of UXOR stars with disk inclinations derived from the 
model, as already described above. Despite these dif- 
ficulties, the model of curved puffed up rims by Isella 
et al. (2006) is quite promising and represents the most 
advanced description of the inner disk region so far. 

4.3. Dust grain size in the inner disk region 

Dust properties of the inner protoplanetary disk can 
also be deduced from the models of observed near IR 

1 Monnier et al. (2006) suggest that curved puffed up rim model 
by Isella & Natta (2005) may be more suitable, but notice that 
this model also produces images with a large skewness (Isella et al. 
2006) 
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visibilities. Since in §3.3.2 we argue that NIR visibili- 
ties of high-L objects are not dust related, we limit our 
discussion on dust grain size to low-L YSOs only. 

As already described in §3.1, the inner radius of opti- 
cally thick multigrain disks is controlled by the largest 
grains in the mix. Also, dust growth and dynamics may 
change the opacity structure of the inner disk and form a 
large optically thin zone populated by big (>l/zm) grains 
extending closer to the star than the optically thick part 
of the disk. Results from §3.1 support models of optically 
thick inner disks containing big grains. Smaller grains 
can hide inside the dusty disk, but we can see only big 
grains that populate the disk surface and shield smaller 
grains from the direct stellar heating. Dusty halo mod- 
els of NIR visibilities also favor micron size grains, but 
no larger than ~10/^m (see §3.2). Visibilities are, there- 
fore, clearly suggesting that circumstcllar dust exists in 
the region where submicron grains cannot survive when 
directly exposed to the stellar radiation. 

Presence of big grains in the inner disk has been sug- 
gested before. Grady et al. (1995) modeled variations in 
the UV and visual spectra of UX Ori and obtained the 
best fit with grains >0.15/im, which suggests consider- 
able grain growth in comparison with the ISM dust, van 
Boekel ct al. (2004) obtained spectra from the 1-2 AU 
zone of the immediate surrounding of several low-L Her- 
big Ae/Be stars using mid IR interferometry. These spec- 
tra are dominated by micron size thermally processed 
(crystalline) dust grains, which differs from the outer (2- 
20 AU) region where grains are smaller (^O.l^m) and 
less crystalline. This is in agreement with results based 
on the NIR visibilities, where we expect such a grain size 
gradient due to dust sublimation. Small grains (<0.1/zm) 
around low-L Herbig Ae/Be stars sublimate at >1 AU, 
while micron size grains can survive within 1 AU. These 
big grains populate the inner disk surface and dominate 
detected visibilities and spectra of the inner disk. 

4.4. The difference between high-L and low-L YSOs 

The clear distinction between visibilities of low-L and 
high-L YSOs may be a sign of more extensive dissimi- 
larities between these two classes of objects. The change 
in circumstellar environment detected by visibilities ap- 
pears at ~T0 3 £ Q , which corresponds to the spectral type 
of about B3-B5 during the PMS phase. Here we point 
out some other significant changes in the YSO properties 
that happen around that spectral type. 

Dispersal of circumstellar gas: The mean gas density 
within a radius of 0.8 pc around high-L Herbig stars de- 
creases by almost two orders of magnitude during their 
fast (< 3 x 10 5 yr: Palla & Stahler 1993) evolution to 
the main sequence (Fuente et al. 1998). In contrast, en- 
vironments of low-L stars experience a decrease of less 
than an order of magnitude. The NIR interferometry 
targets high-L objects that have still not dispersed their 
immediate environment. In §3.3.2 we speculate that disk 
accretion rates under such conditions are high enough 
to enable a gaseous accretion disk to dominate the NIR 
emission and visibilities. This differs from low-L objects, 
where disks have time to evolve to the point where ac- 
cretion rates become low enough to make the gaseous 
disk transparent, allowing dust emission to dominate the 
NIR. 



Disk Mass: Disk masses in high-L Herbig stars are sig- 
nificantly lower than in low-L (Fuente et al. 2003). The 
trend is most noticeable in the ratio of disk to stellar 
mass, but it is also significant when the absolute disk 
mass is considered. It could be that these stars evolve 
through the PMS phase too quickly to accrete a massive 
disk before the star disperses its circumstellar environ- 
ment. 

Radiation pressure: Detailed calculations of radiation 
pressure force on dust grains are very complicated and 
uncertain due to our limited knowledge of geometrical 
properties of individual dust grains (Burns et al. 1979). 
Nevertheless, we can exploit some general properties of 
the ratio of radiation pressure force to gravity, 0, and 
derive an estimate for the size of grains repelled by the 
star (0 > 1). The ratio is generally proportional to 
L+/M+. Differences in the shape of stellar spectrum, 
grain refractive indexes, and grain density can scatter the 
value of by an additional order of magnitude (Burns 
et al. 1979). But luminosity differences between low-L 
and high-L YSOs are so large that, in the first approxi- 
mation, a comparison of values in YSOs is dominated 
by A PMS star of L* = 10 3 L Q has ~ 6M 

(Palla & Stahler 1993), which means ^170 times larger 
than in the Solar System. High-L objects have even 
larger L*/M*, hence, even larger 0. The Sun can eject 
dust particles of ^0.2/im in size (Landgraf et al. 1996; 
Burns ct al. 1979), while other particles smaller than 
^10^m in the Solar System have ~ 0.01 — 1. This 
means that we can safely assume that YSOs of 10 3 i Q 
repel all dust grains smaller than ^lO^tm. In general, 
high-L YSOs have this limit at grain sizes larger than 
10/xm because their luminosity is > 1O 3 L0. It is diffi- 
cult to maintain dusty structures around high-L objects 
with such high values of 0. A dusty halo would be very 
efficiently dispersed. Puffed up disk rims cannot main- 
tain their vertical stability and would be efficiently blown 
away. Only dust within the optically thick gaseous disk 
or strongly dragged by gas toward the star can resist 
the radiation pressure force. The situation is completely 
different for low-L objects. A 10 2 L Q Herbig star has 
L+/M* ~ 30 (Palla & Stahler 1993), hence the upper 
grain size limit of dust with > 1 is in the micron range. 
Curiously enough, this is consistent with smaller grains 
not surviving in the surface of the inner disk (see §4.3 
above). 

Variability: The range of photometric variabilities of 
high-L YSOs is much smaller (AV<1.5 mag) than that 
of low-L YSOs (AV<6 mag) (Rodgers 2001; Herbst 
& Shevchenko 1999; Bibo & The 1991; Finkenzcllcr & 
Mundt 1984). This dependence of variability on lumi- 
nosity is not understood. It may be that high luminosi- 
ties are masking underlying accretion variabilities (un- 
like FU Ori stars where the stellar luminosity is much 
smaller than the accretion luminosity) or that dust sub- 
limates so far away from the star that the probability for 
a transient dust obscuration is highly reduced (Grinin 
& Kozlova 2000). Near IR visibilities, however, indicate 
that differences in the geometry of circumstellar mat- 
ter distribution may be the reason. Another interesting 
property of the variability data is the correlation between 
the NIR flux excess and the range of variability ampli- 
tudes. Herbst & Shevchenko (1999) show this for Herbig 
Ae/Be stars by plotting the variability range in V band 
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against the ratio of NIR excess luminosity to the stellar 
luminosity. Skrutskie et al. (1996) show this for T Tau 
stars by plotting variations of intrinsic (K-L) color ex- 
cess against the amplitude of the K-band variability. If 
the puffed-up inner dusty-disk rim model is responsible 
for V-band variability (which we challenge; see §4.1) then 
this correlation indicates temporal variations in the rim's 
height. The halo model, on the other hand, explains 
this correlation by connecting the amount of dust in the 
halo with the dynamics of dusty outflow. A larger NIR 
flux requires a larger optical depth of the halo, which in 
turn is a result of enhanced dust dynamics supplying the 
halo with dust. Side effects of this enhanced dynamics 
are dust obscuration events and NIR emission variabili- 
ties. Since the existence of such a halo is consistent only 
with low-L YSOs (see §3.2), the variability amplitude is 
smaller in high-L YSOs because they lack the halo and 
its dust dynamics. Unfortunately, the current data are 
not sufficient for establishing which theory is correct. 

5. SUMMARY 

We have reviewed all published NIR (H- and K-band) 
visibility data of YSOs and devised a method for their 
model-independent comparison. The method is based on 
scaling the distance of objects and their luminosity out 
of the measured baseline (equation 1). This removes the 
apparent dependence of the object's size on: i) its dis- 
tance and it) radiative transfer scaling due to luminosity. 
Hence, the visibility dependence on scaled baseline de- 
tects inherent differences in the geometry of circumstel- 
lar matter distribution without applying any additional 
theoretical model-dependent assumption. 

The comparison shows a clear distinction between low- 
L YSOs (L* < 10 3 £ Q ) and high-L YSOs (£* > 10 3 L Q ), 
as already suggested by previous studies (Monnier et al. 
2005; Eisner et al. 2004). Low-L visibilities cluster at 
spatial scales ~7 times larger than scales derived from 
the visibility clustering of high-L YSOs. Next, we ana- 
lyze the observed visibility clusters with three types of 
image models. Modeling the whole data cluster instead 
of individual objects reveals or reaffirms some collective 
properties of these objects. 

The first model is the uniform brightness ring, where 
we use dust sublimation as the boundary condition for 
the ring's inner radius. High-L YSOs are inconsistent 
with dust sublimation and appear much too small. The 
model was also not successful in explaining the size of 
T Tau cluster. These stars appear slightly larger than 
model predictions. The model was more successful in 
low-L Herbig Ae/Be stars where the visibility cluster can 
be modeled with optically thick rings of 0-60° inclination 
and a dust sublimation temperature of ~1500K. The sec- 
ond model is the optically thin dusty halo. It explains 
the T Tau cluster with halos of ~0. 2-1.0 visual optical 
depth. Low-L Herbig Ae/Be stars require optical depths 
of ^0.15-0.8. Halos made of micron size grains provide 
the best fit to these data clusters. Finally, the third 
model is a classical accretion disk. We use this model on 
high-L Herbig Be stars and show that it can accommo- 
date observed visibilities. The model does not constrain 
the accretion rate because equally successful fits can be 
built with = - 10" 4 M Q /yr. We argue, however, 
that the NIR emission from accretion disks must be due 
to gas and not dust, hence, accretion rates must be high 




Fig. 12. — The emerging picture of the inner protoplanetary disk 
structure based on observed properties of high-L and low-L YSOs. 
Top: Low-L objects have two competing models explaining their 
NIR visibilities and anomalously high NIR excess: (A) puffed-up 
inner disk rim and (B) dusty outflow creating a halo around the 
inner disk. Detected variability due to dust obscuration events 
suggests that either (C) the height of puffed-up rim is variable and 
temporarily blocks the view toward the star or (D) clumps of dust 
appear in the dusty outflow and occasionally intercept the line of 
sight. The disk inside the zone of dust sublimation (E) is optically 
thin, while the rest is optically thick (G) due to dust. Gaseous 
stellar and disk wind (F) are also present. Bottom: High-L YSOs 
have a simpler structure. The star is surrounded by optically thick 
gaseous accretion disk (A), which extends much closer to the star 
than the dust sublimation distance. But dust still may survive 
within the optically thick disk interior (B). Intense gaseous stellar 
and disk wind (C), combined with the stellar radiation pressure, 
are efficiently dispersing the surrounding environment. (S) marks 
the star. 



in order to produce required thermal emission from the 
gas. 

We also discuss variability properties of YSOs and 
made an attempt to incorporate them into the exist- 
ing models of inner protoplanetary disk. We argue that 
dust obscuration events detected through photometric 
and spectroscopic variability are caused by dust clumps 
too high above the dense inner protoplanetary disk to be 
considered dynamically part of the disk. Instead, we ad- 
vocate a model where dust is ejected out of the disk and 
blown away by the radiation pressure. This would create 
a clumpy dusty halo-like outflow above the inner disk. 
Unfortunately, we are not aware of any physical process 
that could lift the dust out of the disk. Without this pro- 
cess the only alternative is to assume that the inner disk 
rim is puffed-up and undergoes occasional height varia- 
tions responsible for dust obscuration events. But this 
model is also problematic because it requires disk inch- 
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nations lager than actually observed in objects with dust 
obscuration events. All these facets of evidently compli- 
cated circumstellar environment of YSOs are sketched in 
Figure 12. Grain size also plays a very important role in 
this environment. We argued for large (>l/xm) grains in 
the surface of inner disk region of low-L YSOs based on 
theoretical and observational evidence. 

In addition to the difference in their visibility clusters, 
high-L and low-L YSOs differ in many other observa- 
tional aspects. High-L Hcrbig Be stars: i) are much more 
efficient in dispersal of their circumstellar gas and dust, 
ii) have smaller disk masses, Hi) have strong radiation 
pressure capable of destroying circumstellar dust struc- 
tures made of grains up to lO/xm or more in size, iv) 
are far less photometrically variable. All these proper- 
ties clearly indicate that the geometry of circumstellar 



gas and dust distribution around high-L YSOs is funda- 
mentally different from those in their lower luminosity 
counterparts. Our current understanding of these envi- 
ronments is far from satisfactory and it is largely incom- 
plete. 



We thank Doron Chclouche and Anatloy Mirosh- 
nichenko for useful discussions. We also thank the anony- 
mous referee for useful comments that have improved this 
paper. D.V. acknowledges supports by the NSF grant 
PHY-0503584 and the W.M. Keck Foundation. T.J. ac- 
knowledges the hospitality and financial support of the 
Institute for Advanced Study during his visit to the in- 
stitute. 



REFERENCES 



Acke, B., van den Ancker, M. E., & Dullemond, C. P. 2005, A&A, 
436, 209 

Acke, B., & van den Ancker, M. E. 2004, A&A, 426, 151 
Akeson, R. L., et al. 2005a, ApJ, 635, 1173 
Akeson, R. L., et al. 2005b, ApJ, 622, 440 

Akeson, R. L., Ciardi, D. R., van Belle, G. T., & Creech-Eakman, 

M. J. 2002, ApJ, 566, 1124 
Akeson, R. L., Ciardi, D. R., van Belle, G. T., Creech-Eakman, 

M. J., & Lada, E. A. 2000, ApJ, 543, 313 
Bell, K. R., & Lin, D. N. C. 1994, ApJ, 427, 987 
Berrilli, F., et al. 1992, ApJ, 398, 254 

Beust, H., Karmann, C, & Lagrange, A.-M. 2001, A&A, 366, 945 

Bibo, E. A., & The, P. S. 1991, A&AS, 89, 319 

Burns, J. A., Lamy, P. L., & Soter, S. 1979, Icarus, 40, 1 

Calvet, N., Muzerolle, J., Briccno, C, Hernandez, J., Hartmann, 

L., Saucedo, J. L., & Gordon, K. D. 2004, AJ, 128, 1294 
Calvet, N., Hartmann, L., & Strom, S. E. 2000, Protostars and 

Planets IV, 377 

Calvet, N., Hartmann, L., & Kenyon, S. J. 1993, ApJ, 402, 623 
Cesaroni, R., Galli, D., Lodato, G., Walmsley, C. M., & Zhang, Q. 

2006, Protostars and Planets V, in press (astro-ph/0603093) 
Chavarria-K., C, de Lara, E., Finkenzeller, U., Mendoza, E. E., & 

Ocegueda, J. 1988, A&A, 197, 151 
Chen, C. H., & Jura, M. 2003, ApJ, 591, 267 
Chini, R. 1981, A&A, 99, 346 

Cieza, L. A., Kessler-Silacci, J. E., Jaffe, D. T., Harvey, P. M., & 

Evans, N. J., II 2005, ApJ, 635, 422 
Claria, J. J. 1974, A&A, 37, 229 

Clarke, C. J., & Pringle, J. E. 2006, MNRAS, 370, L10 

Cohen, M., Bieging, J. H., Welch, W. J., & Dreher, J. W. 1985, 

ApJ, 292, 249 
Colavita, M., et al. 2003, ApJ, 592, L83 

Danchi, W. C, Tuthill, P. G., & Monnier, J. D. 2001, ApJ, 562, 
440 

de Lara, E., Chavarria-K., C, & Lopez-Molina, G. 1991, A&A, 
243, 139 

de Winter, D., & van den Ancker, M. E. 1997, A&AS, 121, 275 
Dorschner, J., Begemann, B., Henning, T., Jaeger, C, & Mutschke, 

H. 1995, A&A, 300, 503 
Doucet, C, Pantin, E., Lagagc, P. O., & Dullemond, C. P. 2006, 

A&Aaccepted (arXiv:astro-ph/0608615) 
Drew, J. E., et al. 1997, MNRAS, 286, 538 
Dullemond, C. P., & Dominik, C. 2005, A&A, 434, 971 
Dullemond, C. P., Dominik, C, & Natta, A. 2001, ApJ, 560, 957 
Dullemond, C. P., van den Ancker, M. E., Acke, B., & van Bockcl, 

R. 2003, ApJ, 594, L47 
Finkenzeller, U., & Mundt, R. 1984, A&AS, 55, 109 
Eiroa, C, et al. 2002, A&A, 384, 1038 

Eisner, J. A., Hillenbrand, L. A., White, R. J., Akeson, R. L., & 

Sargent, A. I. 2005, ApJ, 623, 952 
Eisner, J. A., Lane, B. F., Hillenbrand, L. A., Akeson, R. L., & 

Sargent, A. I. 2004, ApJ, 613, 1049 
Eisner, J. A., Lane, B. F., Akeson, R. L., Hillenbrand, L. A., & 

Sargent, A. I. 2003, ApJ, 588, 360 
Friedemann, C, et al. 1993, A&A, 277, 184 
Friedjung, M. 1985, A&A, 146, 366 



Fuente, A., Rodriguez-Franco, A., Testi, L., Natta, A., Bachiller, 

R., & Neri, R. 2003, ApJ, 598, L39 
Fuente, A., Martin-Pintado, J., Bachiller, R., Neri, R., & Palla, F. 

1998, A&A, 334, 253 

Garcia Lopez, R., Natta, A., Testi, L., & Habart, E. 2006, A&A, 

accepted (arXiv: astro-ph /0609032) 
Grady, C. A., et al. 2000a, ApJ, 544, 895 

Grady, C. A., Sitko, M. L., Russell, R. W., Lynch, D. K., Hanner, 

M. S., Perez, M. R., Bjorkman, K. S., & de Winter, D. 2000b, 

Protostars and Planets IV, 613 
Grady, C. A., Perez, M. R., The, P. S., Grinin, V. P., de Winter, 

D., Johnson, S. B., & Talavera, A. 1995, A&A, 302, 472 
Green, J. D. et al. 2006, ApJ, accepted (astro-ph/0605365) 
Grinin, V. P. 2000, ASP Conf. Ser. 219: Disks, Planetesimals, and 

Planets, 219, 216 
Grinin, V. P., & Kozlova, O. V. 2000, Astrophysics, 43, 239 
Grinin, V. P., Kozlova, O. V., Natta, A., Ilyin, I., Tuominen, I., 

Rostopchina, A. N., & Shakhovskoy, D. N. 2001, A&A, 379, 482 
Grinin, V. P., The, P. S., de Winter, D., Giampapa, M., 

Rostopchina, A. N., Tambovtseva, L. V., & van den Ancker, 

M. E. 1994, A&A, 292, 165 
Hartmann, L., & Kenyon, S. J. 1996, ARA&A, 34, 207 
Hartmann, L., Kenyon, S. J., & Calvet, N. 1993, ApJ, 407, 219 
Herbig, G. H., Andrews, S. M., & Dahm, S. E. 2004, AJ, 128, 1233 
Herbst, W., & Shevchenko, V. S. 1999, AJ, 118, 1043 
Herbst, W., Warner, J. W., Miller, D. P., & Herzog, A. 1982, AJ, 

87, 98 

Hernandez, J., Calvet, N., Briccno, C, Hartmann, L., & Bcrlind, 

P. 2004, AJ, 127, 1682 
Hillenbrand, L. A., Strom, S. E., Vrba, F. J., & Kecnc, J. 1992, 

ApJ, 397, 613 

Hofmann, K.-H., Balega, Y., Ikhsanov, N. R., Miroshnichcnko, 

A. S., & Weigelt, G. 2002, A&A, 395, 891 
Hutchinson, M. G., Albinson, J. S., Barrett, P., Davies, J. K., 

Evans, A., Goldsmith, M. J., & Maddison, R. C. 1994, A&A, 

285, 883 

Isella, A., & Natta, A. 2005, A&A, 438, 899 

Isella, A., Testi, L., & Natta, A. 2006, A&A, 451, 951 

Ivezic, Z., Ncnkova, M. & Elitzur, M. 1999, User Manual for 

DUSTY, Univ. of Kentucky internal report, accessible at 

http://www.pa.uky. edu/^ mo she/ dusty 
Ivezic, Z., & Elitzur, M. 1997, MNRAS, 287, 799 
Ivezic, Z., & Elitzur, M. 1996, MNRAS, 279, 1011 
Johns-Krull, C. M., & Valenti, J. A. 2001, ApJ, 561, 1060 
Kenyon, S. J., Yi, I., & Hartmann, L. 1996, ApJ, 462, 439 
Lada, C. J. 1985, ARA&A, 23, 267 

Landgraf, M., Augustsson, K., Griin, E., & Gustafson, B. A. S. 

1999, Science, 286, 2319 
Levreault, R. M. 1988, ApJ, 330, 897 
Liu, M. C, et al. 1996, ApJ, 461, 334 

Lorcnzetti, D., Giannini, T., Nisini, B., Benedettini, M., Elia, D., 

Campeggio, L., & Strafella, F. 2002, A&A, 395, 637 
Lynden-Bell, D., & Pringle, J. E. 1974, MNRAS, 168, 603 
Malbet, F. et al. 2006, A&A, in press 
Malbct, F., et al. 2005, A&A, 437, 627 
Malbet, F., et al. 1998, ApJ, 507, L149 



20 



Vinkovic & .Turkic 



Mannings, V., & Sargent, A. I. 1997, ApJ, 490, 792 

Manoj, P., Maheswar, G., & Bhatt, H. C. 2002, MNRAS, 334, 419 

Meeus, G., Waters, L. B. F. M., Bouwman, J., van den Anckcr, 

M. E., Waelkens, C, & Malfait, K. 2001, A&A, 365, 476 
Millan-Gabet, R., et al. 2006a, ApJ, 645, L77 
Millan-Gabet, R., et al. 2006b, ApJ, 641, 547 
Millan-Gabet, R., Schloerb, F. P., & Traub, W. A. 2001, ApJ, 546, 

358 

Millan-Gabet, R., Schloerb, F. P., Traub, W. A., Malbet, F., 

Berger, J. P., & Bregman, J. D. 1999, ApJ, 513, L131 
Miroshnichenko, A., & Corporon, P. 1999, A&A, 349, 126 
Monnicr, J. D. 2003, Reports of Progress in Physics, 66, 789 
Monnier, J. D. et al. 2005, ApJ624, 832 
Monnier, J. D., et al. 2006, ApJ, 647, 444 
Monnier, J. D. & Millan-Gabet, R. 2002, ApJ579, 694 
Mora, A., ct al. 2004, A&A, 419, 225 

Muzerolle, J., Calvet, N., Hartmann, L., & DAlessio, P. 2003, ApJ, 
597, L149 

Muzerolle, J., D'Alessio, P., Calvet, N., & Hartmann, L. 2004, ApJ, 
617, 406 

Muzerolle, J., Luhman, K. L., Briceno, C, Hartmann, L., & Calvet, 

N. 2005, ApJ, 625, 906 
Natta, A., Grinin, V. P., Mannings, V., & Ungerechts, H. 1997, 

ApJ, 491, 885 

Natta, A., Grinin, V. P., & Tambovtseva, L. V. 2000, ApJ, 542, 
421 

Natta, A., & Whitney, B. A. 2000, A&A, 364, 633 

Natta, A., Testi, L., & Randich, S. 2006, A&A, 452, 245 

Palla, F., & Stahler, S. W. 2000, ESA SP-445: Star Formation from 

the Small to the Large Scale, 179 
Palla, F., & Stahler, S. W. 1993, ApJ, 418, 414 
Perez, M. R., van den Anckcr, M. E., de Winter, D., & Bopp, B. W. 

2004, A&A, 416, 647 
Pogodin, M. A., Malanushenko, V. P., Kozlova, O. V., Tarasova, 

T. N., & Franco, G. A. P. 2006, A&A, 452, 551 
Porter, J. M., Drew, J. E., & Lumsden, S. L. 1998, A&A, 332, 999 
Prusti, T., & Mitskevich, A. S. 1994, ASP Conf. Ser. 62: The 

Nature and Evolutionary Status of Herbig Ae/Be Stars, 62, 257 
Przygodda, F., van Boekel, R., Abraham, P., Melnikov, S. Y., 

Waters, L. B. F. M., & Leinert, C. 2003, A&A, 412, L43 
Quanz, S. P., Henning, T., Bouwman, J., Ratzka, T., & Leinert, 

C. 2006, ApJ, accepted (astro-ph/0605382) 
Rodgers, B. M. 2001, Ph.D. Thesis, 

Rodgcrs, B., Wooden, D. H., Grinin, V., Shakhovsky, D., & Natta, 

A. 2002, ApJ, 564, 405 
Rostopchina, A. N., Grinin, V. P., & Shakhovskoi, D. N. 2001, 

Astronomy Reports, 45, 51 
Sandell, G., & Weintraub, D. A. 2001, ApJS, 134, 115 
Shakhovskoj, D. N., Grinin, V. P., & Rostopchina, A. N. 2005, 

Astrophysics, 48, 135 
Shepherd, D. S., & Churchwcll, E. 1996, ApJ, 472, 225 



Shevchenko, V. S., Ibragimov, M. A., & Chenysheva, T. L. 1991, 

Soviet Astronomy, 35, 229 
Sitko, M. L., Carpenter, W. J., Lynch, D. K., Russell, R. W., 

Grady, C. A., Brafford, S. M., & Wooden, D. H. 2004, American 

Astronomical Society Meeting Abstracts, 205, #16.06 
Skrutskie, M. F., Meyer, M. R., Whalen, D., & Hamilton, C. 1996, 

AJ, 112, 2168 
Straizys, V., et al. 1989, A&A, 222, 82 

Strom, K. M., Strom, S. E., Breger, M., Brooke, A. L., Yost, J., 
Grasdalen, G., & Carrasco, L. 1972, ApJ, 173, L65 

Strom, S. E., Grasdalen, G. L., & Strom, K. M. 1974, ApJ, 191, 
111 

Tuthill, P. G., Monnicr, J. D., Danchi, W. C, Hale, D. D. S., & 

Townes, C. H. 2002, ApJ, 577, 826 
Tuthill, P. G., Monnier, J. D., & Danchi, W. C. 2001, Nature, 409, 

1012 

Terranegra, L., Chavarria-K., C, Diaz, S., & Gonzalez-Patino, D. 

1994, A&AS, 104, 557 
Tcsti, L., Natta, A., Shepherd, D. S., & Wilner, D. J. 2003, A&A, 

403, 323 

Testi, L., Palla, F., & Natta, A. 1998, A&AS, 133, 81 

van Boekel, R., Waters, L. B. F. M., Dominik, C, Bouwman, J., 

de Koter, A., Dullemond, C. P., & Paresce, F. 2003, A&A, 400, 

L21 

van Boekel, R., et al. 2004, Nature, 432, 479 

van Boekel, R., Min, M., Waters, L. B. F. M., de Koter, A., 
Dominik, C, van den Ancker, M. E., & Bouwman, J. 2005, A&A, 
437, 189 

van den Anckcr, M. E., Blondel, P. F. C, Tjin A Djie, H. R. E., 

Grankin, K. N., Ezhkova, O. V., Shevchenko, V. S., Guenthcr, 

E., & Acke, B. 2004, MNRAS, 349, 1516 
van den Anckcr, M. E., Wesselius, P. R., & Tielens, A. G. G. M. 

2000, A&A, 355, 194 
van den Anckcr, M. E., de Winter, D., & Tjin A Djie, H. R. E. 

1998, A&A, 330, 145 
Vinkovic, D. 2006, ApJ, 651, 906 

Vinkovic, D., Ivezic, Z., Miroshnichenko, A. S., & Elitzur, M. 2003, 

MNRAS, 346, 1151 
Vinkovic, D., Ivezic, Z., Jurkic, T., & Elitzur, M. 2006b, ApJ, 636, 

348 

Yudin, R. V. 2000, A&AS, 144, 285 

Warren, W. H., & Hesser, J. E. 1978, ApJS, 36, 497 

Weigelt, G., Balega, Y. Y., Hofmann, K.-H., & Preibisch, T. 2002, 

A&A, 392, 937 
White, R. J., & Ghez, A. M. 2001, ApJ, 556, 265 
Wichmann, R., Bastian, U., Krautter, J., Jankovics, I., & Rucinski, 

S. M. 1998, MNRAS, 301, L39 
Wilkin, F. P., & Akeson, R. L. 2003, Ap&SS, 286, 145 



